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ABSTRACT 



Context. 

Aims. The optical variability of the blazar PKS 2155-304 is investigated to characterise the red noise behaviour at largely different 
time scales from 20 days to 0(> lOyrs). 

Methods. The long-term optical light curve of PKS 2155-304 is assembled from archival data as well as from so-far unpublished 
observations mostly carried out with the ROTSE-III and the ASAS robotic telescopes. A forward folding technique is used to de- 
termine the best-fit parameters for a model of a power law with a break in the power spectral density function (PSD). The best-fit 
parameters are estimated using a maximum-likelihood method with simulated light curves in conjunction with the Lomb Scargle 
Periodogram and the first-order Stmcture Function (SF). In addition, a new approach based upon the so-called Multiple Fragments 
Variance Function (MFVF) is introduced and compared to the other methods. Simulated light curves have been used to confirm the 
reliability of these methods as well as to estimate the uncertainties of the best-fit parameters. 

Results. The light curve is consistent with the assumed broken power-law PSD. All three methods agree within the estimated uncer- 
tainties with the MFVF providing the most accurate results. The red-noise behaviour of the PSD in frequency / follows a power law 
with f~P, p = 1 -8^2 and a break towards f° at frequencies lower than f m 
Conclusions. 



(2.7^ yrs)" 1 . 
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1. Introduction 

It is widely accepted that the broad-band emission of an ac- 
tive nucleus of a galaxy (AGN) is largely dominated by non- 
thermal (synchrotron and inverse Compton) emission from rel - 
ativistic electrons or positrons (or protons) (lUlrich et al.lll997h . 
In the optical band, thermal emission from the central region 
(broad or narrow-line regions, accretion disk) and stellar emis- 
sion from the host galaxy additionally contribute. In blazars (to 
which PKS 2155-304 belongs) the jet emission is commonly 
considered to be dominant over all other (thermal) contributions 
consistent with the absence of prominent emission or absorption 
lines in the optical spectrum. While it is still an open question 
how and if changes in the accretion rate propagate into varia- 
tions of the jet emission, there are two common explanations for 
the variabilit y to be observed from the jet: self-organised criti- 
cality (SOC) (iBak et al J 1987.). either from the accretion process 
or jet process resulting in a re d noise behavi our, and jet dynamics 
causing periodic oscillations dRiegerll2004l) . 

Over sufficiently long time scales the average flow of mat- 
ter drifting from the outer region of the accretion disk to- 
wards the inner most stable orbit is believed to be constant. On 
smaller time scales, the matter may fall i n avalanches of dif - 
ferent size and duration towards the centre (iDendv et al.lll999l) . 
iKawaguchi & Mineshigd (1 19991) showed that this can be simu- 
lated in a simple automaton model producing SOC. If the opti- 
cal emission correlates to the released energy of the avalanches, 



this model results in a red noise light curve. A more general mo- 
tivation on SOC in BL Lac objects is given by the models of 
SOC in compact pla smas (ISivronll9 98) and astronomical shocks 
dMalkov et al.ll2000t) . At long time scales the effect of SOC may 
be limited by the effective size of the emitting region. Therefore, 
it is expected that at longer time scales the SOC and therefore 
the red noise behaviour vanishes. 

PKS 2155-304 is an optical bright X-ray sel ected BL Lac 
type o bject (my a 14 ma s) at a redshift of z = 0. 1 16 dFalomo et al.l 
Il993l) . Ground based observations in the near-infrared have been 
used to separate the relatively faint ( M r = -24.4 mag ) ellipti - 
cal host galaxy fro m the bright nucleus jKqtilainen et al.ll 1998b . 
ISmith et all (1 19921) and iHeidt et al.l (1 19971) observed a varia- 
tion of the optical flux of 2% and 10% respectively per day. 
iPaltani et al.l (1 19971) found fast variations on the time scales of 
f m in < 15 min which were limited by the noise of the instruments. 
They also found that on a time scale of days the variability is 
compatible with red noise with a power-law index of /3 — 2.4. 
X-ray observations indicate a similar power-law behaviour rang- 
ing from 1 . 7 to 3.5 during se veral continuous pointings lasting 
up to a day dGaur et al. 1 120101) . Recently, red noise has been de- 
tected from PKS 2155-304 at very high energies (> 100 GeV) 
on timescales of hours u p to frequencies of ^ H z in the power 
spectral density function dAharonian et al . 2001. 

Searches for periodiciti es in the optical revealed evidences 
of a periodicity of 0.7 days (lUrrv et al.lll993l) . 
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In this paper a detailed « 7 years long light curve of 
PKS 2155-304 measured with the ROTSE-III telescope sys- 
tem is presented. Archival and almost simultaneous observa- 
tions are added to the data to obtain a long-term light curve 
covering timescales up to as 77 years. Subsequently the intrin- 
sic power spectral density function (PSD) is analysed with the 
Lomb Scargle Periodogram, the Structure Function and a new 
method called Multiple Fragments Variance Function. The re- 
sult is compared to the results of simulated light curves with red 
noise behaviour with different shapes of the PSDs. The best-fit 
parameter set for the assumed model is found with a maximum- 
likelihood analysis. 

In Section [2] a summary of the observations is given. In 
Section [3] the methods and maximum-likelihood analysis used 
to characterise the long-term behaviour of the optical flux of 
PKS 2155-304 are described in detail. Readers who are mainly 
interested in the results, may go directly to Section|H where the 
long-term light curve of the object and the results of the meth- 
ods and maximum-likelihood analysis are presented. The paper 
closes with a conclusion in Section[5] 



2. Summary of observations 

2.1. The ROTSE telescope system and data analysis 

The ROTSE-III telescope systerrQ consists of four robotic tele- 
scopes (called ROTSE-IIIa, b, c and d) located in Australia, 
USA, Namibia, and Turkey. The telescopes are of a Ritchey- 
Chretien type with a parabolic primary mirror (® 450 mm) and a 
focal length of 850 mm. The field of view encompasses 1.85° x 
1.85°. Only light with a wavelength between 400 and 900 nm 
can pass the broad-band filter. On average each telescope col- 
lects 250 observational frames per night. Typically, the limit- 
ing magnitude is tss 17.3 mag (up to 19 mag at perfect conditions) 
and a few thousand individual objects are detected within each 
frame. ROTSE-III is equipped with a fully automated computer- 
controlled data acquisition syst em (DAQ) as well as an au- 
tomatic data-reduction pipeline. lAkerlof et al.l d2003l) provides 
more details on the system. 

The observations used for this analysis have been carried out 
whenever PKS 2155-304 was observable. Observations are taken 
as pairs of exposures of 20 or 60 s each with a dithering pointing 
pattern to reduce the impact of hot pixels. Pairs of observations 
are taken at latency intervals of approximately 45 minutes. 

Each frame is dark-field subtracted and flat-field corrected. 
The thin entrance window makes an additional fringe-pattern 
correction necessary. After the objects in the frame have been 
extracted by means of the SExtractor (lBertinll2002h software, 
their coordinates and instrumental (uncalibrated) magnitudes are 
matched and adjusted to the USNO A2.C0 catalogue. During 
this process, individual objects are flagged as problematic (e.g. 
blended point-spread functions, saturated pixels etc.). 
The photometric data of all objects detected in a sequence 
of frames taken on a field are combined in so-called match- 
structures. During this process, a relative correction of fluctua- 
tions in the photometry is applied, which benefits from the large 
field of view with a correspondingly larger number of detected 
objects. De tails on the data calibr ation and analyses chain are 
provided in lRvkoff & Smith! d2003l) . 

Starting August 24th, 2003, PKS 2155-304 has been ob- 
served regularly with the ROTSE-IIIc telescope system in 



http : //www . ro tse.net/| 

http : // tdc-www . harvard . edu/catalogs/ua2 . html 



Table 1: Data log of the ROTSE observations. 

Total number of frames 12 946 

Number of frames with sufficient quality 6 310 
Average flux in mjy 
Standard deviation in mJy 



15.9 
6.9 



Table 2: Overview of the archival data used in this work. 



Period 



Number of data points Reference 



1934 - 1940 

1979 - 1995 

1980 - 2007 
1998 - 1999 



61 
134 
1721 
12 



(1) 
(2)... (13) 
(14)... (16) 

(17) 



Refer ences. 

( fl979l). (3) iMiller & McAlister! 



AIL 



Hamuv&Mazd j 1987), 
Mead et aU 



IGriffiths et alj dl979l) : 

' ( 119831) . (4 ) iBrindle et al.1 



(2) 



IGriffiths et al.1 
' (119861). 



(6) iPica et al . 1 1988) . (7) | Treves et al.1 



(1990), (9) 



(lO)ICarini & Milled dl992T>. ( ll )ISmith etaLldl992l) . (12) 
d!993l> and (13) lUrrv et al.1 d!993 |) as collected in 
Il996l): (14) ICampbelll d2004h . (15) ICourvoisier eta! 
il997h 



Smith &Sitkd d!991l). 



Jannuzi et al 



(16) lPesceetal.1 d!997 

(17) lTommasi et alJd200ll) . 



Zhang & Xie 
'I995fr and 



collected" in lOsterman et al* I d2007l) 



Namibia. Additional observations have been performed with 
ROTSE-IIIa in Australia since October 10th, 2003. Data up to 
October 26th, 2010 are presented and analysed in this work. 
Table Q] shows statistics of the data on PKS 2155-304. To as- 
sure high data quality and reduce systematic errors, images with 
inferior quality reported by the DAQ are rejected. The complete 
light curve is presented in Section l4~Tl Fig. |4] A table containing 
the data is available at the CDS. Column 1 lists the Julian date 
of the observation, Column 2 and 3 give the apparent luminosity 
(R-Band) with the statistical uncertainty. 

Since ROTSE does not take use of an R-band filter but cal- 
ibrates the instrumental magnitude to the catalogue values of 
stars, the measured R-magnitudes of PKS 2155-304 suffer from 
a systematic deviation due to the non-star like spectrum. This 
requires a cross calibration with simultaneous R-band observa- 
tions taken with other instruments (e.g. the ATOM telescope in 
lAharonian et al . 200§. 

2.2. Archival data and observations 

Archival (mostly photometric plate) data on PKS 2155-304 cov- 
ering the period from 1934 to 2007 have been collected from 
several publications as listed in Table [2] 

Furthermore, 397 mostly simultaneous observations have 
been obtained w ith the All Sky Automated Survey (ASAS) 
(|Poimanskill2002l) from November 25th, 2000 till October 28th, 
2008. ASAS is a project dedicated to continuous photometric 
monitoring of approximately 10 7 stars brighter than 14 magni- 
tude to search for transient and variable objects. The data used 
here were taken with the ASAS-2 instrument located at Las 
Campanas, Chile. The ASAS-2 consists of a / = 135 mm, 
f/D = 1/1 .8 lens system and a Pictor 416 with a KODAK KAF- 
0400 CCD sensor (512x768 pixel). ASAS has been monitoring 
the southern hemisphere since 2000. 

Where required the magnitudes given for a different wave- 
length band are adjusted to the R-band by adding a constant 
offset obtained from observa tions taken with the ATOM tele- 
scope ( lAharonian et al.ll2009l) . The apparent luminosities of the 
whole data set (ROTSE, ASAS and archival data) are con- 
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verted into flux units with 13.3 mag corresponding to 10 mJy 
(lAharonian et alJl2005h . 



2.3. Statistical and systematic uncertainties of the 
photometry 

The estimated statistical error is crucial for simulating the de- 
tector noise behaviour, and therefore we have carried out a care- 
ful analysis of the detector noise as a function of the observed 
brightness. For the ROTSE-III data, we have estimated the total 
(statistical and systematic error) from the data: based upon 660 
individual pointings on PKS 2155-304, we have selected w 3 000 
other objects that are covered in these pointings. For each of 
these 3 000 light curves, the square root of the variance o~ rotse 
has been calculated. The correlation between ctrotse and the 
flux <J> of an object follows a power-law shape which has been 
parametrised with 



o~ rotseC^) 
mJy 



0.0764 



— ) 

lOmJy/ 



1.302 



+ 0.0564. 



(1) 



In order to neglect any variable objects median values have been 
used to perform a least-squares fit. The median values are calcu- 
lated from the values within bins of 1 mJy. The values and the 
fitted power-law are shown in FigJT] 




10 20 30 

Average flux in mJy 



40 



Fig. 1: Standard deviations ctrotse of the light curves of » 3 000 
individual objects represented by the data points. The red curve 
is the least-squares fit to the medians of the values within 1 mJy 
bins. 



For the ASAS data, the quoted photometric uncertainty has 
been parametrised as a function of the total flux: 



mJy 



= 0.162 



lOmJy 



1.128 



(2) 



The photometric uncertainty on the archival plate data has 
been quoted to be in the range of 0.01 ... 0.1 magnitudes. These 
values are used to simulate the absolute uncertainties for the cor- 
responding simulated data points. 



3. Time series analysis 

Two well-established methods - the Lomb-Scargle Periodogram 
(LSP) (see Appendix lA.il ) and the Structure Function (SF) (see 
Appendix I A. 21 ) - as well as a newly developed method called 
Multiple Fragments Variance Function (MFVF) (see Section 
13.11 ) are used to characterise the timing behaviour of the optical 
long-term light curve of PKS 2155-304. 



3.1. Multiple Fragments Variance Function (MFVF) 



Start. 



Split light curve into two 
fragments at the centre. 



Treat fragment as light 
curve. 



Count number of data points 
in the second fragment. 




Calculate the variance and 
the length of the fragment. 



Count number of data points 
in the first fragment. 




Omit this fragment. 



yes 



Omit this fragment. 



Fig. 2: Flowchart of the Multiple Fragments Variance Function 
(MFVF). 



We have developed an alternative method to characterise the 
variability of the source at different timescales in a simple and 
robust way. The variance 



var(a) := (((a) 



a?) 



(3) 



of the entire light curve is computed and subsequently, the light 
curve is subdivided at the half of the time lengths of the interval 
into two fragments and the variance for each of the two indepen- 
dent fragments is computecB This process of subdividing the 
light curve continues recursively until the number of remaining 
measurements within the fragment drops below a critical value 
(see Fig.[2]for a flowchart of this process). For the analysis car- 
ried out here, the requirement on the minimum number of mea- 
surements within the fragment is chosen to be > 10. This number 
is a compromise between the limited number of measurements 
and the smallest timescale resolvable. 

The fragment length is calculated by the time difference be- 
tween the first and last measurement available within this frag- 
ment. The variance is therefore calculated for light curve seg- 
ments which cover intervals with widely different values of dura- 
tion. Finally, these values are averaged (binned) within 100 bins 



3 Alternatively the interval may be split at the position of the largest 
gap- 
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per decade of fragment length. This function yields information 
about the variance of the variability depending on the time scale. 

Like the SF the MFVF is computed in the time domain. 
Unlike the other methods it is not sensitive to periodicities. 
It rather concentrates on information about the mean variabil- 
ity amplitude of the light curve occurring on a specific time 
scale. Therefore it is a goal-oriented method to find a maxi- 
mum time scale of the variability. The sampling of the light 
curve is taken directly into account, so the distorting effect of 
uneven sampling is limited to an acceptable level (see below). 
An IDL as well as a Python script for the MFVF can be found at 
https : //sourcef orge . net/pro j ects/mfvf / 

3.2. Comparison of the methods with simulated light curves 

Since all methods suffer at various degrees from the effect of un- 
evenly sampled light-curves (sampling effects), we use sets of 
simulated light-curves, in order to take sampling effects and the 
influence of instrumental noise into account. These simulated 
light-curves are useful to compare the behaviour of the three 
methods. 

The simulated light curves are generated under the assump- 
tion of a red noise with a power spectral density function (PSD) 
following a power-law (PSD °c f~P above a characteristic fre- 
quency / m in). The original recip e for generating pure red-n oise 
light-curves first introduced by iTimmer & Koeni j d 19951) has 
been extended to include a break at f m \ n . For frequencies below 
fmin, we assume a white noise behaviour {PSD = const.). The 
normalisation of the simulated light curve is chosen to match the 
mean flux level and the variance of the measured light curve. 
In order to investigate sampling effects the light curves were 
generated either with an even sampling (ideal case) or with an 
uneven sampling (realistic case). For both cases the total length 
A?total °f the light curve and the number of data points N \, s are 
chosen to match the values of the light curve of PKS 2155-304. 
For the ideal case the observation times are spread evenly over 
Motai- This results in an evenly sampled light curve with one 
measurement every Af to tai/M>bs = 3.37 days and without any 
gaps. 

For the realistic case the sampling is chosen to be exactly 
the same as it is for the light curve of PKS 2155-304. The raw 
evenly sampled simulated light curves are sampled down to the 
actual measurement times. Only for the realistic case each indi- 
vidual flux measurement has been randomised with a Gaussian 
function in order to account for experimental uncertainties on 
the measurements (see Section 12.31 for details on these uncer- 
tainties). For each data point the variance of the Gaussian errors 
were chosen individually to simulate the according instrumental 
noise. 

For each case (ideal and realistic), 5 000 independent light 
curves (covering a 27 000 days each) have been simulated with 
ft = 2.0 and / rain = 10~ 3 d" 1 . For the ideal case, all methods 
should produce results which are close to the intrinsic behaviour 
of the light curve. Here, the sampling effects are limited to the 
periodic sampling and to the finite length of the simulated light 
curves. By comparing the results from the realistic and ideal case 
simulations, sampling effects related to the uneven sampling can 
be investigated. 

The individual values of the methods are accumulated in 
two-dimensional histograms. In Fig. [3] the histograms for the 
LSP (upper panel), SF (middle panel), and MFVF (lower panel) 
for the ideal case (left side) and the corresponding histograms 
for the realistic case (right side) are displayed. The resulting 



probability-density function (PDF) for a fixed value of frequency 
is shown in logarithmic bins (10 per decade). 

The ideal case LSP (Fig. [3] upper panel left side) reproduces 
the simulated shape of a broken power-law with a visible devi- 
ation at high frequencies (> 0.1 days -1 ) related to the sampling 
frequency of 0.297 days 1 . The LSPs from the uneven sampled 
light curves (right side) show a considerable flattening deviat- 
ing from the original red noise behaviour. Additionally, the po- 
sition of the break is smeared out over roughly half a decade 
in frequency. The flat (white) noise behaviour is distorted as 
well, showing a decrease of the LSP for decreasing frequencies. 
Besides the broad band distortion of the LSP, features in nar- 
row frequency bands appear as a consequence of the sampling 
effects. Most notably is the drop in the LSP at a frequency of 
1 year 1 due to the quasi-periodic sampling. The ideal case SF 
(Fig. [3] middle panel, left side) follows the expected behaviour 
for red noise ST t" with a « ft - 1 dHughes et al. I [1991 
lLainela & Valtaojalll993l) . For values of t > 0.3/ mm the SF flat- 
tens out indicating the break in the simulated power-law. The 
break in the SF is smeared out into a gradual roll-over. The real- 
istic case does not deviate systematically from the ideal case. 

Finally, the result of analysing the light curve using the 
MFVF shows a consistent behaviour between the ideal and re- 
alistic simulations. In both cases, the simulated break at f . = 

' J min 

1 000 d translates into a bend in the slope of the MFVF at a frag- 
ment length of w 1 000 d. However for the ideal case no frag- 
ments shorter than « 40 d are found where the number of obser- 
vations is greater than 10 due to the even sampling. Therefore 
the MFVF is not defined at those time scales. 

3.3. Best-fit estimates for / min and ft 

Given that the sampling of the light-curve is generally not even, 
the results of the methods described above suffer from sampling 
and instrumental effects. Therefore, it is necessary to apply a 
forward folding technique to reliably reconstruct / mm and ft. 
We have simulated 5 000 realistic light curves (as described in 
Section [3T2b for each parameter set on a grid of values of 0.8 < 
ft < 7.2 in steps of Aft = 0.2 and -4.8 < log 10 (/ min /d -1 ) < -2.2 
increasing in steps of 0.2. For each set of parameters we have 
generated the PDFs for the LSP, SF and MFVF as described in 
Section [3~!2l For the measured light curve a likelihood-estimator 

N 

mfnm) = Xlog/Mft/ata), (4) 
(=1 

is calculated with p,(ft, /min) representing the probability to mea- 
sure a value of e.g. the £SP(fi) for a given set of parameters 
ft, /min in the i-th bin. p : {ft,fm in ) is looked-up directly in the 
corresponding PDF(/?, / m i n ). The best-fit PDF(ft, / m ; n ) maximises 
■Cift, /min) and determines the best-fit parameter set. It is note- 
worthy that the PDFs are not necessarily symmetric which re- 
quires the use of a likelihood in stead of a least-squares function 
dEmmanoulopoulos et alj|2010l) . However, given that the indi- 
vidual probabilities are not independent of each other, the like- 
lihood function is not well-defined. Nevertheless, the parame- 
ters which maximise the likelihood are good estimators for the 
best-fit parameters. It is not straightforward to estimate the un- 
certainties on the best-fit parameters given that the shape of the 
likelihood function (but not the position of its maximum) are 
considerably modified by the correlations between the individ- 
ual probabilities (see also Section l4~3l on tests of the method and 
estimates for the uncertainties of the parameters). Since all sim- 
ulated light curves are normalised to the same mean flux level 
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10 



10" 3 10" 2 10 

Frequency / in days" 1 



Time lag t in days 
10 3 10 2 10' 




10" 3 10" 2 10" 

Frequency / in days" 1 
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10' 
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Fragment length t in days 
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10"" 10" 3 10" 2 10" 

Frequency / in days" 1 
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Fig. 3: Two dimensional histograms of LSP, SF, and MFVF (from top to bottom) values of artificial light curves with fi = 2.0 and 
/min = 10~ 3 d . On the left side histograms are shown for the case of an even sampling of one measurement every 3.37 days while 
on the right side, the effects of a realistic (uneven) sampling and simulated instrumental noise are visible (see text for details). 
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and variance of the measured light curve, the methods are not 
sensitive to these parameters. 



4. Results 

4.1. Long-term optical light curve of PKS 2155-304 

The optical R-band light curve measured with ROTSE is shown 
in Figure|4] The contribution of the host galaxy (15.1 mag ) is sub- 
tracted here and in the following. The gaps in the light curve are 
caused by bad weather periods or when PKS 2 155-304's position 
is too close to the sun. The light curve consists of 6 310 indi- 
vidual data points. The peak to peak amplitude observed during 
this time interval varies from 5 to 40 mJy indicating a relative 
variability of 35 mJy or a factor for the flux of « 8. The variabil- 
ity amplitude is largest when comparing flux states separated by 
long time intervals (few years) while it decreases when compar- 
ing flux states on shorter timescales. This behaviour is typical of 
red noise. 

In addition to the ROTSE-III measurements archival and 
contemporaneous data on PKS 2155-304, as described in 
Section [2] have been included in Fig. [5] Overall, the light 
curve covers the time from 1934 to 2010 with an extended gap 
between 1940 and 1979. Apparently, the maximum amplitude of 
variability in the long-term light curve is similar to the observed 
variability in the seven years of ROTSE-III observations. Before 
analysing the light curve in detail, this is a first indication that 
the red noise behaviour does not extend towards time scales 
much longer than a few decades. 



4.2. Results of the time series analysis 

We have applied the three different methods described in 
Sections fA.lllA.2H3.ll to analyse the variability of the long-term 
light curve shown in Fig. [5] Based upon the resulting LSP, SF, 
and MFVF, the best-fit estimators for the power-law index (3 as 
well as the break frequency f m \ n have been derived individually 
following a log-likelihood method described in Section l3~3l 

The instrumental noise for the archival data might be under- 
or overestimated. Therefore the variability on small time scales 
of a simulated light curve could be distorted to an unrealistic 
case. So we use only bins at time scales > 20 days for the log- 
likelihood method to minimise systematic errors arising from 
this uncertainty. 

The LSP of the long-term light curve of PKS 2155-304 is 
plotted in the top panel of Fig. [6] together with the PDF (his- 
togram in colour scale) of the LSPs of the simulated light curves 
with the best-fitting parameters. The shape of the measured LSPs 
is in good agreement with the PDFs derived from simulations. 
As expected, sampling effects lead to a considerable distortion 
of the intrinsic power-law shape. Additionally, the uneven sam- 
pling produces features in the LSP which again are reproduced 
in the simulated data-sets. 

In a similar fashion, the SF (middle panel in Fig-Eb as we H as 
the MFVF (bottom panel in Fig- |6j are shown together with the 
PDF for the best fit. For both cases, the simulated and measured 
data-sets agree well with each other. Both, broad band shape and 
narrowband features from sampling effects are well reproduced 
in the simulated data. 

For all methods, the resulting best-fit parameters for /3 and 
/min with their lcr uncertainties are listed in Table [3] The uncer- 
tainties for the parameters are estimated as described in Section 



Table 3: Best-fit parameters estimated by the methods separately. 



Method 


P 


logcwd-') 


LSP 
SF 
MFVF 


5 +n 

1-6-85 


r> o+O.l 
J-A-0.4 

_3 4+0.4 

-3 o +0 - 3 



14.31 Graphs with the log-likelihood functions for the parameters 
can be found in AppendixlBl 

4.3. Reliability of the method and estimates of the 
uncertainties 

The likelihood method introduced above suffers from complex 
correlations between individual probabilities. This in turn re- 
quires to confirm that the method of maximising the likelihood 
estimator does actually provide a meaningful estimator of the 
underlying true parameters. Furthermore, it is necessary to esti- 
mate the statistical uncertainties on the best-fit parameters with- 
out relying on the usual method of calculating the second deriva- 
tives of the likelihood in the maximum. 

Therefore an empiric approach is indicated. We have tested 
the reliability of the likelihood estimation with simulated light 
curves created with certain (well known) parameter sets. Again, 
5 000 realistic light curves (as described in Section I3.2l i have 
been simulated with the best-fit parameter set which has been 
found with the respective method. The light curves are analysed 
with the same methods as applied to the real light curve result- 
ing in a distribution of best-fit parameter sets. The resulting dis- 
tributions are then used to verify that the simulated parameters 
are reconstructed reliably and the width of the distributions pro- 
vide a robust estimate of the uncertainties on the parameters. In 
Fig.|7]the two-dimensional histograms are displayed as contours 
(from top to bottom LSP, SF, and MFVF). The cross marks the 
simulated (input) parameters and the error bars indicate the 1 cr 
confidence region for each parameter of the marginalised distri- 
bution^ which are plotted along the axes. 
For the SF and MFVF, the input parameters are well recovered 
with relatively small uncertainties (see Table [3]l while for the 
LSP, the reconstructed values are spread out over a large interval 
of f3 resulting in large errors. While / m ; n is in very good agree- 
ment with the results of the SF and the MFVF, /? differs notice- 
ably. However, this is in accordance to the large uncertainties of 
/3 for the LSP, which is obviously affected by sampling effects. 
The LSP's uncertainties for /3 in the direction for higher values 
is probably underestimated, because the border of the grid is too 
close to the maximum of the likelihood. 

It remains to be demonstrated that the model and its best-fit 
parameters are a good characterisation of the data. Again, we 
use an approach relying on 5 000 simulated light curves to check 
the consistency of data and assumed model. As a goodness of fit 
estimate the value of £, of the observed light curve is compared 
to the distribution of £ of the light curves simulated with the 
best-fitting values of j3 and / m j n . For the three different methods 
a fraction of 4 % (LSP), 3 % (SF), and 50 % (MFVF) of the 
simulated data are fit worse by the model than the actual data. 
Therefore the model is a reasonable description of the variability 
of the optical flux of PKS 2155-304. 

4 The 1 cr intervals have been calculated by integrating over one of 
the two parameters. The resulting one dimensional histogram is then 
normalised and the interval is truncated by taking 15.85% from both 
sides symmetrically such that the central interval contains 68.3 % of the 
distribution 
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Fig. 4: Optical light curve (R-band) of PKS 2155-304 measured with ROTSE. 
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Fig. 5: Optical long-term light curve (R-band) of the core region of PKS 2155-304. The light curve combines ROTSE observa- 
tions as well as other a rchival data (solid bla ck circles: ROTSE, op en purple circles: A ll-Sky-Automated-Sys tem, red open squares: 
Osterman et al.l (120071) . green solid squares: IZhang & Xid (119961) . blue solid circles: [Griffiths et al.l d 1979b . black open squares: 
Tommasi et aTT(l200lh V 



5. Conclusion 

We have presented an optical light-curve covering seven years 
of densely sampled (6 310 individual pointings) observations of 
PKS 2155-304 with the ROTSE-III telescopes in Namibia and 
Australia. Combined with archival data and ASAS observations, 
a long-term (sparsely sampled) light-curve going back to 1934 
is analysed to characterise the red noise behaviour under the as- 
sumption that the power spectral density function (PSD) follows 
a power-law with slope (3 and a break at a minimum frequency 
/min- At frequencies below / m j n , a white noise (/? = 0) PSD is 
assumed. For the analysis only time scales > 20 days have been 
considered to avoid systematic errors due to uncertain estimates 
of the instrumental noise for the simulated light curves. The 
three different methods for analysing the long-term variability 
of PKS 2155-304 produce consistent results for simulated light 



curves with /? « 1.8 and f m [ n « 1/1 000 d _1 . The most precise 
result (discussed in the following) is given by the newly intro- 
duced Multiple Fragments Variance Function (MFVF), which is 
a straightforward method in the time domain to reveal informa- 
tion about the mean variability amplitude depending on differ- 
ent time scales. The result obtained with the first-order Structure 
Function (SF) is in very good agreement with the result of the 
MFVF. The slope of the Lomb-Scargle Periodogram (LSP) is 
affected strongly by sampling effects, so the uncertainties of the 
best-fit parameter [3 are considerably larger for this method than 
for the MFVF and SF. 

The power-law in dex of the PSD P = 1-8^2 i s harder than 
the value obtained by iPaltani et al.l (1 19971) at timescales smaller 
than 10 d (/? = 2.4+02). This may be an indication for an ad- 
ditional break or curvature in the PSD on short timescales. The 
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power-law index of the red noise observed at other energy scales 
are similar. This hints at a common origin of the observed emis- 
sion at the base of the jet. 

Our analysis confirms the presence of red noise up to 
timescales of a few years without significant indications for any 
additional features as e.g. quasi-periodic oscillations. The exis- 
tence of a break in the PSD at / m ; n = 1/1 000 d -1 is found to be 
significant at the level of w 6 cr judging from the error on / m j n es- 
timated for the MFVF method. This break corresponds to a time 
scale of T max * 2.7*j'? yrs after which the red noise behaviour 
disappears. In the context of SOC this time scale is linked to 
the "memory" of the system. At time scales longer than r max the 
properties of the emission process de-correlate from its past. Due 
to causal reasons, the spatial dimension of the emission region is 
constrained to be r < cr max = 0.84*9 - 5? pc. 

The observed red-noise from PKS 2155-304 extends over 
five orders of magnitude in the time domain and over twelve or- 
ders of magnitude in the photon energy with apparently similar 
power-law slope. This implies that there is probably a link con- 
necting particle acceleration processes on spatial scales ranging 
from O(100 AU) to 0(1 pc). 
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Fig. 6: Normalised Lomb Scargle Periodogram (LSP), Structure 
Function (SF) and Multiple Fragments Variance Function 
(MFVF) (from top to bottom). The green lines/points are the 
measured results of the light curve of PKS 2155-304. The his- 
tograms in colour scale represent the PDFs obtained from simu- 
lated light curves with the best-fit parameter sets (J3 = 5.0, / m ; n = 
10- 3 - 2 d- 1 ), 08 = 1.6, / min = 10- 3 - 4 d- 1 ) and 03 = 1.8, / mi „ = 
10 3 d ') respectively. Only the bins at / < 0.05 d 1 are used 
for the maximum-likelihood method. 
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Fig. 7: Distribution of maximum likelihood estimates of sim- 
ulated light curves with parameters found for PKS 2155-304. 
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are calculated to obtain the 1 cr error bars. 



9 



M. A. Kastendieck et al.: Long-term optical variability of PKS 2155-304 



Appendix A: Standard methods of time series 
analysis 

A. 1. Lomb-Scargle Periodogram (LSP) 

The normalised LSP (lLombll 1 9761: IScarddfl98l is a method to 
estimate the true underlying PSD of the variability of unevenly 
sampled time series data. For a given set of values a-, := a(ti) 
measured at the times t> the LSP is defined as 



ZSP{u) = 

2cr- 



1 f \Lj(aj-a) cos(aj(tj - e))] 2 



{ Z,cos 2 (w(/;-e)) 

[Zj(fly-g) sin (^( f ;-e))]" 
2,sin 2 (a>(tj-e)) 



(A.l) 



with a> = luf, cr 2 being the variance and a the mean value 
of the time-series a,. The time delay 



1 _, Zjsinh.cotj) 

— tan - 

2w cos (2u>tj) 



(A.2) 



is required to achieve that £S'P(a)) is invariant under a change 
of the phase dPress & Rvbickilll989l) . In order to weigh widely 
different frequency scales equally, we calculate the LSP for val- 
ues off with a constant interval of 0.01 decades on a logarithmic 
scale. 



A.2. Structure Function (SF) 

The first order SF dSimonetti et al .11 1985b for a time-lag r be- 
tween any pair of observations is defined as 

ST{t) = - V [a & + r)-a fc)] 2 , (A3) 
n 

i, 

with n being the number of summands. The structure function 
can only be computed for discrete values of r that exist in the 
time series: 



t e f,- - tj\ : ti,ti e T\ 



(A.4) 



where T is the set of times where measurements have been car- 
ried out. Commonly, the values of the structure function are 
combined and averaged in intervals of r to smooth fluctuations 
on time scales shorter than the width of the bins. In order to 
weigh widely different time scales equally, we choose a constant 
bin width of 0.01 decades on a logarithmic scale. 

Appendix B: Log likelihood 

The log-likelihood values for the grid of the parameters are 
shown in the panels (from top to bottom LSP, SF, and MFVF) 
of Fig. EH] The log-likelihood functions £ of the SF and MFVF 
have a similar shape with a ridge-like region with a bend. For the 
SF and the MFVF the maximum is located in the region where 
the ridge bends. For the LSP, the maximum of the likelihood is 
located on a ridge which is oriented along the direction of in- 
creasing p, which obviously does not constrain the best-fit value 
of very well. 
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Fig. B.l: Logarithm of the likelihood for the real light curve be- 
ing realised by red noise with a certain fi and f m \ n . Top for the 
LSP, middle for the SF and bottom for the MFVF. The peaks 
located in the white shaded areas are marked with a black cross. 
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